A standard Hamiltonian formulation for the dynamical Casimir effect 
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We present a quantum description of piioton creation via dynamical Casimir effect based on tiie 
standard Hamiltonian formulation. The particle representation is constructed in the expansion of 
field operators fixed with the initial modes. The Hamiltonian is presented in terms of the creation and 
annihilation operators with the time-varying couplings which originate from the external properties 
such as an oscillating boundary or a plasma mirror of a semiconductor slab. Some consideration is 
also made for the experimental realization with a semiconductor plasma mirror. 
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Introduction.-The quantum nature of vacuum provides 
a variety of physically interesting phenomena, including 
the Casimir effect [l]. The so-called dynamical (non- 
stationary) Casimir effect (DCE), as well as the static 
force, has been investigated extensively B, H, 3, H, @, 0, 
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references therein), where photons are created from the 
vacuum fluctuation in non-adiabatic change of the sys- 
tem driven by vibration of a cavity or expansion of the 
universe. Most of the theoretical approaches are based on 
the field expansion in terms of the instantaneous modes. 
Since there is no unitary map among the instantaneous 
modes with different boundary conditions, the Hamilto- 
nian formulation to describe the time-evolution does not 
exist in the standard sense Even in this case, the 
time-evolution of the instantaneous-mode operators may 
be described by means of an effective Hamiltonian 
On the other hand, a suitable transformation may be 
made to move to a specific coordinate system with fixed 
boundaries for canonical quantization [sl, [^, 0, [22| . 

Experimentally, it is difficult to realize a sufficient mag- 
nitude of mechanical vibration at a resonant frequency 
~ GHz to create a significant number of photons for de- 
tection. As a feasible alternative, it has been proposed 
recently that the oscillating wall can be simulated by a 
plasma mirror of a semiconductor slab irradiated by pe- 
riodic laser pulses 15 1 . (See also Refs. [H, [l^].) 

In this paper, we investigate a quantum description 
of DCE, presenting the standard Hamiltonian to gov- 
ern the unitary time-evolution. We are particularly con- 
cerned with the experimental realization of DCE with a 
semiconductor plasma mirror. It is indeed important to 
develop the Hamiltonian formulation to investigate the 
quantum properties of the system, including the detec- 
tion of created photons through interaction with suitable 
probe such as atoms. This standard description has sev- 
eral advantages: (1) The particle representation is con- 
structed in the expansion of field operators fixed with 
the initial modes. It is neither necessary to trace the 
mode change in time, nor to seek a specific coordinate 
system for quantization. Simple formulas are presented 
to calculate the couplings among the creation and anni- 



hilation operators in the Hamiltonian as space-integrals 
involving the initial mode functions. Similar formulas are 
obtained for the instantaneous-mode effective Hamilto- 
nian where the time-derivatives of the mode functions are 
further involved (2) The time -variation of the cre- 

ation and annihilation operators in the standard formula- 
tion precisely represents the unitary quantum evolution 
all the time. On the other hand, the time-variation of 
the instantaneous-mode operators and that of the mode 
functions together provide the quantum evolution. The 
instantaneous-mode operators and mode functions coin- 
cide with those of the standard formulation just at each 
period of the oscillation. (3) The present formulation is 
applicable to various physical setups, including the oscil- 
lating wall and the semiconductor plasma mirror, as in- 
vestigated later. We can check that under the situations 
where the mode functions do not change largely in time, 
as usually considered, this standard description provides 
essentially the same result for DCE as the instantaneous- 
mode description. Furthermore, the Hamiltonian is read- 
ily calculated even for the large time-variation of exter- 
nal properties, clarifying the dependence on experimental 
parameters specifically for the plasma mirror case. 

Standard Hamiltonian formulation.-We consider a 
scalar field generally in 3+1 space-time dimensions. The 
Lagrangian is given as 
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(1) 



-(V</>)2--m2(x,t)^2 

(n = c = 1) 0, i, 0, [ii, [li. Here, e(x, t) and m^{:ic, t) 
represent the dielectric permittivity and conductivity (ef- 
fective "mass" term), respectively, in the matter region 
such as a semiconductor slab. As specified later, they are 
space-time dependent, simulating the boundary oscilla- 
tion. Conventionally, the instantaneous modes (x, t) 
(real, orthonormal and complete) at each time t with 
time- varying frequencies uJa(t) are adopted according to 
the boundary oscillation: [— -I- m^(x, t)]/Q,(x, t) = 
e(x,i)w2(i)/^(x,t) with e(x, t)/„(x, t)/^(x, i)^^^ = 
5af}/\^u}a{t)]- Instead, we here construct the particle rep- 
resentation in terms of the initial modes 



/^(x) = /„(x,i = 0), u:l=uj^{t = Q). 



(2) 
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The canonical field operators in the Heisenberg picture 
are expanded with the creation and annihilation opera- 
tors at,(i) and aa{t) as 

0(x,t) = Y.M)+^imi{^), (3) 

a 

n(x,t) = e(x,0)^^c.°[-a„(^)+4(^)]/°(x), (4) 

a 

where n(x, t) = dC/d(j) = e(x, t)(/)(x, i). Then, the 
Hamiltonian (Schrodinger picture) is presented by the 
usual procedure as 



Hit) 



1 f 

y2le(x,t) 



+ <I>[-V^ + m^{x,t)](t)\d^x 



a,/3 



(5) 



where the space-integral is taken over the whole region V 
which is "fixed" (not time-dependent) according to the 
physical setup, as explicitly shown later for typical cases. 
The Heisenberg and Schrodinger pictures are related by 
the unitary transformation U{t) of time-evolution gener- 
ated by this Hamiltonian H{t) as aa{t) = W{t)aaU{t). 
Hence, the time-dependence of the couplings in Eq. (O, 
which originates from the c-number external quantities 
e(x, t) and m^(x, t), is common in any pictures related by 
unitary transformations. The Hamiltonian in the Heisen- 
berg picture Hult) — U'^ {t)H{t)U{t) is obtained simply 
by ^ aa{t), a^a ^ ^alO; 4> ^ 0(x,<), H ^ H(x, i) in 
H{t). 

The mode frequencies ijJa{t), intermode couplings 
/ic(^(i) and squeezing terms gapit) are calculated by con- 
sidering the orthonormality of (x) which obey the wave 
equation with e(x, 0) and to'^(x, 0): 



iOa(,t) 

gap{t) 



-*[-G^^(t)+G^'^(t)], 



(6) 
(7) 
(8) 



SV{t) 



sv{t) eA(x, i) 
mi(x,i)/0(x)/0(x)d3a: 



/0(x)/^"(x)d3^, (9) 
(10) 



The integrals for G^'™(t) are evaluated in practice in the 
subregion 5V{t) (C V), possibly time-dependent, where 
e(x, t) and m^(x, t) vary in time as e^"'^(x, t) = e~^(x, t) — 
e-i(x,0)andmi(x,t) EETO2(x,i)-m2(x,0). [G^'"(0) = 
at t = with H{Q) diagonalized in terms of /°(x).] In 
order to demonstrate the relevance of this Hamiltonian 
formulation for DCE, we investigate two typical cases in 
the effective 1-f 1 dimensions (1) oscillating wall and (2) 
plasma mirror of a semiconductor slab. 



Oscillating wall.-The boundary walls may be repre- 
sented by high potential barriers of matter extending in- 
finitely (or finite and long) outside the cavity as 

m^(x,t) = m^[-oo < X < 5{t),L < x < oo]. (11) 

Here, the right side is fixed at x — L, while the left side 
varies in time around x = as 6(0) = < S{t) < Si ^ L. 
The dielectric is taken uniformly as e{x, t) — ei in the 
matter region. This setup with potential barriers, rather 
than the rigid boundary conditions, may be similar to 
Ref. [1^ where the matter-field interaction is treated 
dynamically. The mode functions are given as 



/fc(x,t) 



(j^\k'\(x-5(t)) 

Asinfc[a; — 5{t) 

QQ-\k'\{x-L) 



(— OO, 5{t)) : wall 

-i] m,L] (12) 

(L, +oo) : wall 



with the dispersion relations = (fc^ + fc^)/eo = (fc'^ + 
fc^ -|-m^)/ei, where k'^ ^ — < and m ~ \k'\ ^ k ^ 
1/L for the large m^, and kj_ (^ k) is the momentum in 
the orthogonal spatial 2 dimensions [l2, 13, 21 1. 

For the large potential barrier m^, the frequency mod- 
ulation Su!k{t) and the diagonal squeezing coupling gkk{t) 
in the instantaneous- mode Hamiltonian [9|, |lO| are calcu- 
lated by noting sinfc[i — S{t) + ^] ~ at a; = L with 
^ ~ 1/m < L and \B/A\, \G/A\ ~ fc^ < 1 as 

fcfc(t) ~ iol[S{t)/L]rk, gkkit) ^ di;k{t)/[iLdk{t)], (13) 

where = fc^/eo(a^^)^. (The dielectric contribution is 
suppressed significantly by eiU)\/w? <^ 1.) This leading 
result is independent of the large . By taking formally 
the limit m — > cx) (^ — > 0), the usual moving boundary 
conditions fk{5{t),t) = fk{L,t) = are reproduced. 

On the other hand, in the present standard description, 
the frequency modulation and the squeezing term are 
calculated in Eqs. ©-(Uni) with eo(w°)^/m2 < 1 as 



SuJk{t) ^ 2igkk{t) ^ 



m 



sin^ k(x + £)dx 



dujk{t) 



[m6{t) < 1] 



SLUk{t)[md{t)]^/3 [md{t) > 1] 



(14) 



where (x) 

s{o) = 0, e 



Asmk{x -|- ^) in < X < L with 
1/m, A ~ {Leo(^k)~^^^ (normalization). 



and ■m'ji^{x,t) = m? in 5V{t) = (0, ^(a;)). In the limit 
m ^ oo the standard 5ujk{t) oc [m5{t)Y diverges except 
at i = with 5(0) = 0, while the instantaneous-mode 
fcfc(t) remains finite as Eq. (|13p . This corresponds to the 
claim that the Hamiltonian does not exist in the moving 
boundary problem . The squeezing couplings gkk (t) in 
Eq. (fTS]) and gkk{t) in Eq. (fT4|) appear to be different 
even with SaJkit) ~ dujk{t). It, however, will be shown 
that they provide essentially the same result for the pho- 
ton creation at the resonance. 
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Plasma mirror.-We next consider the case of plasma 
mirror which is realized with a semiconductor slab irra- 
diated by periodic laser pulses [l5[- The dielectric re- 



sponse of plasma is given by e{uj) = ei[l — {ujp/uj'^)] with 
the plasma frequency ujp — (nee^/eim,)^/^ in terms of 
the effective electron mass m* and the conduction elec- 
tron number density rie proportional to the laser power 
Wiaser- The dispcrsion relation in plasma k'^ = e(a;)w^ = 
eiw^ — (riee^/m*) can be taken into account in the slab 
region [I, I + S] around x = I with a thickness (5(<C L) as 

e{x, t) — ei(i), m^{x, t) = rn^{t) = ne{t)e'^ /m^, (15) 

where ^^(0) — for Wiascr(O) = 0. (The spatial distribu- 
tion of the conduction electrons may also be considered 
readily.) The mode functions are given as 



fk{x,t) 



Dsmkx [0,Z) 

Q^^k'x ^_ c-e-^fe'^ [/, ? -I- (5] : slab 

AsiTLk[x - 5 + S,{t)] {l + 5,L\ 



(16) 



{k' = i\k'\ for fc'^ < with large m^). The Dirichlet 
boundary condition (corresponding to the TE mode) is 
adopted at x — Q,L with sinfc[L — 5 + ^{t)] = 0. 

The standard 5LjJk{t) and gkk{t) are calculated in Eqs. 
dSD-dlOl) with Eq. dill) for f{[x) at t = as 



gkk{t) = (^/2)^2[ 



'5ra{t)]/L, 

Mt) + 5,n{t)yL. 



(17) 

(18) 



Here, the effective wall oscillation is enhanced as 



SMI 5 
5,n{t)/6 



-[ei(0)/eo][l - ei(0)/ei(t)] sin^ kl, (19) 



[ml{t)/e,{u;lf]sinHl. 



(20) 



square 



oc 



This effect is almost proportional to the 
of mode function around the slab [/fe(0]^ 
sin^fc; since Sl'^^ [fl{x)fdx 2± [ft{l)fS for k'5 - 
[ei(0)/eo]i/2((5/i) < 1 at < = 0. If the slab is placed 
at the boundary x = I = Q, sir? kl is replaced with 
(fc(5)^/3 ^ {S/LY <^ 1, as observed in Ref. [l9| claiming 
that DCE is suppressed in the TE mode. The signifi- 
cant photon creation, however, will take place even in 
the TE mode if the slab is placed apart from the bound- 
aries X — Q,L which are the nodes of /fc(a:) [Til. [23l|. 

The shift £^{t) in the instantaneous modes of Eq. 
(fTHll is determined mainly proportional to 5 to give the 
frequency modulation 5u)k(t). The squeezing coupling 
9kk{t) is then calculated with the formulas for the ef- 
fective Hamiltonian [13] ■ After some calculations we 
find again the relations 5u)k{t) — SuJk{t) and gkk{t) — 
[i/2iDk{t)]gkk{t), as seen in Eqs. (fT3)) and p4)) for the os- 
cillating wall, where the change of dielectric is assumed 
to be small, \ei{t) ~ ei(0)| < ei(0), as usual This 
ensures the same result on DCE in both the descriptions 
for the case of small oscillation, as shown later. 



The above calculations of duk{t) and gkkif) are valid 
even for the large ei(i) and m^lt) to provide the enhanced 
displacement |(5e^m(i)| 3> 5, which will be plausible exper- 
imentally. It is not necessary here to consider the large 
deformation of the mode functions in time which inval- 
idates the usual perturbative calculation assuming the 
small change of the instantaneous modes. 

Photon creation as squeezing. -Owce the Hamiltonian 
is presented in terms of the creation and annihilation 
operators, the quantum properties of the system are in- 
vestigated readily by using the methods of quantum op- 
tics. We here consider the quantum evolution for DCE, 
restricted to a single resonant mode with time-varying 
frequency ijj{t) = wq + 5uj{t) and squeezing coupling g(t), 
omitting the mode index "fc". The intermode couplings 
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will not provide significant contributions 
since generally due to the non-equidistant frequency dif- 
ferences they are highly oscillating in the rotating-wave 
frame (interaction picture) where the term (aj)a^a is elim- 
inated for the average frequency {u)) = wq + over the 
period T = 2tt/VI of the laser pulse. 

The Heisenberg equation ia{t) — [a{t) ^ H'ii{t)] is de- 
scribed as the master equation, 

A = -iuj{t)A + 2g{t)B., B = iuj{t)B + 2g*{t)A, (21) 

in terms of the Bogoliubov transformation, 

a(t) = A{t)a + B*{t)a\a){t) = A*{t)a) + B(t)a. (22) 

The solution is expressed as A{t) ~ coshr(t)e''^'*'^*\ 
B{t) = sinhr(i)e"^«(*), ensuring \A{t)\^ - \B{t)\'^ = 1 
with A{0) = 1,5(0) = 0. The unitary time-evolution is 
then given as a phase rotation and squeezing, 

U{t) — e«-f(*)g-[^*(t)'"i-A(t)ata+]/2gj0A(t)a+Q ^23) 

withA(t) =r(t)e'['^'^W-'^«(*)l ;2l. The phase factor e*-^(*) 
with K{t) = (f)A{t)+ /p uj{t')/2dt' is included to reproduce 
the zero-point energy of H{t) in iU{t) = H{t)U{t). 

An analytic solution for A{t) and B{t) is obtained in 
the rotating-wave approximation by replacing uj{t) — > 
Wo + (average), g{t) (5)ne~*^* (Fourier compo- 
nent). By noting the time-evolution of the number oper- 
ator a'^{t)a{t) = \B{t)\'^aa'^ -I- . . ., we obtain the photon 
creation via DCE (vacuum squeezing) as 

N^{t) = (0|at(t)a(t)|0) ~ (|2(g}n|/x)' sinh^ (24) 

with the effective squeezing 



X = V|2(5)o|2-A2, (25) 

allowing for the detuning A of the laser pulse [l3, [13] as 

9. = 2(^0 + {5lo) + A). (26) 

The resonance for DCE is given precisely by = 
2{ujo + {5^)) rather than = 2a;o, as considered in the 
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instantaneous-mode approach [l8|. If = 2ujo is taken 
naively with A = —{duj), the effective squeezing x is sig- 
nificantly reduced, even possibly becomes imaginary with 
Nj{t) < 1 oscillating as sin^ \x\t. 

We have solved numerically the master equation typi- 
cally with du!{t) = {Su!){l — cosfit) and g{t) — ~iSuj{t)/2 
to confirm that the rotating-wave approximation is fairly 
good for |(5a;(i)| ^ ojq. The instantaneous- mode solu- 
tion is also obtained with Soj{t) — Suj{t) and g{t) = 
[i/2uj{t)]g{t), as seen so far. It almost reproduces the 
rotating-wave approximation, smoothing the actual small 
oscillation of Nj{t) due to that of (5a; (t). The rela- 
tions Soj{t) = Su!{t) and g{t) = [i/2uj{t)]g{t) really im- 
ply 1 2(g) a I ~ |2(g)n| (Fourier components) around the 
resonance = 2{ujo + {Suj)) for |(5(jj(t)| ^ ujq, giving es- 



sentially the same Nj{t) in Eqs. 



and 



We now discuss the experimental realization of DCE 
with the semiconductor plasma mirror. It will be fea- 
sible with the sufficient maximal laser power W^i"scr 
to achieve the enhanced displacement as S',^^^ ~ 
[(n™'^'^eVeo"^*)/wo]<5 ^ lO^S or larger with sin^ fc^ = 1 
(the slab placed in the middle of cavity I — L/2). In 
this case, the conductivity effect Sm in Eq- ((20| dom- 
inates over the dielectric effect in Eq. (fT9| with 
ei(0) 1 — 10 and ei(0) < |ei(t)| [even for the complex 
ei(t)]. Then, we estimate roughly x(A = 0) = |2(g)o| ~ 
wo(5^''''/i) 10-2cjo for S - lOAtm and L - 0.1m. This 
requires A^puiso ^ 100 repetitions of laser pulse to cre- 
ate iV^ > 10 photons with xi^puiscT) > 1. The cavity 
Q value is reasonable a,s Q > luq/x ^ 10^. The tuning 
of for the resonance should be made with the average 
shift (Suj) ~ ;^ ~ 10~^a;o. The time-profile of Wiaser(*) 
should also be chosen suitably to optimize the Fourier 
component (5)ae~'^* in g{t). A detailed analysis will be 
made elsewhere based on the present formulation. The 
time- varying dielectric function ei(t) (complex) and con- 
ductivity TOp (t) are actually given depending on the laser- 
power profile Wiasor(^)- By using these €i{t) and m'^{t), 
the frequency shift 5uj(t) and squeezing coupling g{t) are 
determined in Eqs. p7|) and (fT8|l . Then, the master 
equation is solved to obtain the photon number Ny(t). 

Detection.-The photons created via DCE can be de- 
tected suitably by Rydberg atoms with principal quan- 
tum number n « 100 and transition frequency ^ GHz 



111 . I23j . Rydberg atoms as two- level system are initially 
prepared in the lower level, and injected into the cavity. 
Some of these atoms are excited to the upper level by 
absorbing the photons, and detected outside the cavity 
as the signal of photons. Recently, high-sensitivity mea- 
surement of blackbody radiation has been performed at 
a frequency 2.527 GHz and low temperatures 67 mK - 
1 K by employing a Rydberg-atom cavity detector with 
a newly developed selective field ionization scheme for 
n « 100 (the atoms excited by absorbing photons are 
selectively ionized by applying an electric field) [23|. It 
exceeds the standard quantum limit, detecting less than 



one photon on average in the cavity. Hence, the single- 
photon detection with Rydgerg atoms is really capable of 
observing even a small number of DCE photons. When 
TVpiyd atoms are injected in the cavity, the number of pho- 
tons detected by atoms is limited roughly as < -/VRyd 
(acutally iVRyd ~ 100 Q). We also note that in order 
to observe purely the vacuum squeezing via DCE, the 
cavity should be cooled well below 100 mK to suppress 
the thermal photons as TV*''^™^! ^ ^ 
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